Resonant microwave properties of a voltage-biased single-Cooper-pair transistor 
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We consider the microwave dynamics and transport properties of a volt age- biased single-Cooper- 
pair transistor. The dynamics is shown to be strongly affected by interference between multi- 
ple microwave-induced inter-level transitions. As a result the magnitude and direction of the dc 
Josephson current are extremely sensitive to small variations of the bias voltage and to changes in 
the frequency of the microwave field. 
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I. INTRODUCTION 

It has long been known that transport properties of a 
superconducting weak link may be strongly affected by 
external time dependent (ac) fields. To give an exam- 
ple, just after the discovery of the Josephson effect it 
was found that an alternating rf-field will rectify the ac 
Josephson current that results from a dc voltage bias and 
produce voltage steps — Shapiro steps — in the current- 
voltage characteristics^. This effect can be viewed as 
a synchronization of the Josephson current oscillations 
of frequency uuj = 2eV/h, generated by a bias volt- 
age V, and the current oscillations induced by the ac 
field. If the Josephson coupling energy Ej is a simple 
harmonic function of the superconducting phase differ- 
ence Ej{<&) = Ej cosQ, the positions of the 
Shapiro steps are given by the expression eV^ 71 ' = nhtu 
(n = 1,2,...), where u> is the frequency of the ac field. 
Deviations from a simple harmonic relation between Ej 
and <f> result in additional "fractional" steps at voltages 
e y(nm) _ n ^_ w ^ TO However their heights are exponen- 
tially small when m n. 

New opportunities to influence the transport proper- 
ties of Josephson junctions by external ac fields appear 
in mesoscopic systems, where only few quantum states 
are responsible for the superconducting charge transfer. 
For such weak links the Josephson energy — no longer 
a macroscopic quantity (being proportional to the area 
of the junction) — may be of the same order of magni- 
tude as the energy quantum Ku of the external ac field. 
If so, the alternating field may significantly affect the 
current by inducing resonant inter-level transitions be- 
tween charge carrying states in the junction and thereby 
changing their population. The internal structure of the 
mesoscopic weak links determines the physical origin of 
the discrete energy spectrum and the relative positions 
of the energy levels inside the superconducting gap 2A. 

Superconducting point contacts (SPCs) and single- 
Cooper-pair transistors (SCPTs) are two fundamentally 
different examples of mesoscopic weak links. In an SPC 
such as the one sketched in Fig. ^l, the rapid variation 
of the superconducting phase across the junction gives 
rise to quasiparticle bound states (Andreev states) that 
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FIG. 1: Schematic energy spectrum of (a) a superconducting 
point contact (SPC) and (b) a single-Cooper-pair transistor 
(SCPT). Energies are plotted as a function of the supercon- 
ducting phase difference $ across the system. An important 
difference between these two examples of a mesoscopic weak 
link is that the discrete energy levels inside the gap A are 
well separated from the continuum levels (dark bands) for an 
SCPT but not for an SPC (see text). 



are localized to the vicinity of the contact. The dis- 
crete energy levels (Andreev levels) associated with these 
states appear inside the superconducting gap as shown 
in Fig. \T£l. Their positions ti arc defined by the single 
electron scattering matrix and in the case of low trans- 
parency contacts one finds that A — |ej| <C A. In an 
SCPT, which is composed of two low-transparency tun- 
nel junctions in series forming a small island as sketched 
in Fig.^D, the Coulomb blockade phenomenon 2 strongly 
affects the structure of the energy spectrum. If the super- 
conducting gap A is greater than the Coulomb energy of 
a grain charged by a single electron, the ground state of 
the SCPT with a completely isolated central dot, for cer- 
tain gate voltages, is degenerate with respect to adding 
one more Cooper pair to the dot^. Weak Cooper-pair 
tunnelling that removes this degeneracy gives rise to cur- 
rent carrying Cooper-pair resonant states, separated — 
as illustrated in Fig.^3 — from the continuum spectrum 
by an energy of order A^iLLS^. 

There is a qualitative difference between an SPC and 
an SCPT in how the energy spectrum depends on the 
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FIG. 2: Sketch of a single-Cooper-pair transistor (SCPT) 
composed of two low-transparency tunnel junctions in series 
forming a small superconducting (S) island bridging the gap 
between two superconducting bulk leads. The potential of 
the island is controlled by a voltage V a applied to the gate 
electrode, which is capacitively coupled to the island. A bias 
voltage V is applied across the SCPT giving rise to a linear 
time dependence of the phase difference «J> = ujjt = 2eVt/h 
between the leads 

superconducting phase change across the system. In an 
SPC the subgap level (as described above this is a discrete 
energy level in the gap between the ground state and the 
continuum excitation spectrum) is mixed with the con- 
tinuum spectrum at certain values of the phase differ- 
ence, while the Cooper-pair resonant levels in an SCPT 
are always separated from the continuum, forming a well 
defined two-level system. The microwave dynamics of a 
biased SPC was explored in Ref. TjJ where it was shown 
that quantum interference among a sequence of tempo- 
rally localized inter-level transitions results in nontriv- 
ial transport properties of a Josephson weak link of this 
type. In an SPC, however, coherence may be preserved 
only during a fixed length of time, since the discrete An- 
dreev states die away as they periodically merge with the 
continuum spectrum at certain values of the phase dif- 
ference. (We recall that a voltage bias makes the phase 
difference increase linearly with time). In contrast, since 
the discrete energy levels are here always separated from 
the continuum, the microwave-induced transport proper- 
ties of a biased SCPT is to a great extent defined by the 
coherent dynamics of the two-level system and strongly 
affected by interference processes. In order to demon- 
strate this phenomenon, we explore below the microwave 
dynamics and transport properties of a biased SCPT. 



II. MODEL 

A sketch of the SCPT structure to be considered here 
is presented in Fig. [21 where a nanoscale superconduct- 



ing island — capacitively coupled to a gate electrode — 
is shown to serve as a bridge between two bulk supercon- 
ducting leads. The Hamiltonian for the system is 

H = ^±91 _ e L j cos($/2 -$)- Ef cos($/2 + fa , 

(!) 

where the Cooper-pair number operator— h and the 
phase operator <f> refer to the superconducting conden- 
sate on the grain, [<j>, ft] — i; Ej 1 "^ is the Josephson 
coupling energy between the grain and the left (right) 
lead, C is the total mutual capacitance and Q — C g V g is 
an induced charge controlled by the gate voltage V g . The 
first term in the Hamiltonian governs the electrostatic en- 
ergy associated with the number of extra Cooper pairs on 
the grain, while the second and third terms describe the 
Josephson coupling between grain and leads. Below we 
consider the symmetric case for which E>j = E$ = Ej. 
Therefore, the Josephson part of the Hamiltonian simpli- 
fies to -2Ej cos(<&/2) cos0. 

It is well known that if the characteristic Coulomb en- 
ergy E q ~ (2e) 2 /2C associated with charge fluctuations 
due to Cooper pair exchange between grain and leads, is 
significantly larger than the Josephson energy Ej , while 
at the same time the Coulomb energy ~ e 2 /2C related to 
single electron fluctuations is smaller than the supercon- 
ducting energy gap, the low-energy states of the grain are 
given by a superposition of only two charge states that 
differ in charge by 2e. Below we will assume that the gate 
voltage is chosen in such a way that Q — — e + SQ where 
\5Q\ <C e. In that case the state of the grain is a superpo- 
sition of the charge state |0), describing the neutral grain, 
and the charge state |1) with one extra Cooper pair. If 
SQ = the Coulomb energy in the states |0) and |1) are 
the same. This degeneracy is, however, removed by the 
Josephson coupling that induces Cooper pair exchange 
between the mesoscopic grain and the bulk electrodes. 
As a result the eigenstates of the Hamiltonian |Q with 
eigenvalues E± — ^fEj cos $ /2 are given by symmetric 
and antisymmetric superpositions of the states |0) and 
|1) : |±) = (|0) ± |1))/V2. These states, which we will re- 
fer to as Cooper-pair resonant states, will — if occupied 
— carry currents given by j± = ±eEjfi~ l sin$/2. 

A bias voltage makes the phase difference over the 
junction change with time, $(t) = uijt + <&o. This time 
dependence generates a periodic variation of the relative 
positions of the Josephson levels as well as a periodic 
change in the current that can be carried by the corre- 
sponding resonant Cooper pair states (see Fig. EJ). The 
total Josephson current through the junction depends 
on the relative "population" Sp = p + — p_ of the lev- 
els and is given by the formula jj = eEj?i~ 5 p sin <f>/ 2 
(p± are probabilities to find the system in the states |±), 
p + + p_ = I). One can find important differences be- 
tween the ac Josephson effect in such a system and in 
ordinary Josephson junctions. If the Coulomb blockade 
is totally lifted, there are no inter-level transitions as the 
levels cross each other and consequently Sp is controlled 
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system, strictly speaking, is not periodic. This is because 
the phases of the complex elements of the scattering ma- 
trix depend on the phase loU + tp which defines the value 
of induced charge at time U. This fact is of great impor- 
tance for the transport properties of the system. As we 
will show below, if the ratio oj/Q, is an irrational num- 
ber the dynamics of the system is qusiperodic and the dc 
Josephson current is equal to zero. But if this ratio is a 
rational number, then a dc Josephson current does flow 
through the junction. The maximum value of the current 
appears when to/Cl is an integer. 



FIG. 3: Time evolution of the Josephson levels E±(t) in the 
single-Cooper-pair transistor (SCPT) of Fig.|3 Dashed lines 
represent microwave- induced inter-level transitions occurring 
at times ti, when the resonance condition \E+(ti) — E-(ti)\ = 
Tiuj is satisfied. 



only by relaxation processes. Even if the temperature is 
much lower than the characteristic inter-level distance, 
a weak relaxation to the ground state with rate v <C luj 
leads to an equalization of the level populations. As a re- 
sult the Josephson current will be proportional to v and 
in the limit v — > not only the dc but also the ac current 
through the Josephson junction vanishes. 

The situation is different if an oscillatory voltage is 
applied to the gate. A time dependent gate voltage 
makes the induced charge vary in time as well, 8Q = 
<5Qo cos(ujt+ip). This in turn gives rise to inter-level tran- 
sitions which control the level populations and by this 
means one may stimulate the Josephson current. The 
most interesting case, which will be considered below, 
appears when the frequency of the gate voltage oscilla- 
tions oj is of the same order as the frequency of charge 
fluctuations on the grain Ej/%. In this case the inter- 
level transitions have a resonant nature and significantly 
influence the Cooper pair transport even for relatively 
small variations of the gate potential. 

In the present work we consider a situation where the 
energy quantum tku of gate voltage oscillations, is of the 
same order but smaller than the maximum inter-level en- 
ergy distance 2Ej/h. At the same time we take u> to be 
much larger than the rate of the gate-voltage induced 
resonant inter-level transitions, e6Qo/hC. Under these 
conditions and if one considers the time evolution of the 
system at low bias voltages, loj <C u>, one can neglect 
inter-level transitions except during short time intervals 
Sti near the times U, when the resonance condition 12 
hui = ±2Ej smflti is fulfilled (£1 — uij/2). During these 
time intervals Landau-Zcner type inter-level transitions 
will occur (see Fig.^J. Therefore the evolution of the sys- 
tem we are interested in is given by a sequence of scat- 
tering events separated by periods of "free" evolution, 
during which the level populations do not change. 

Despite the fact that the level positions vary periodi- 
cally in time with period Tq = 2tt/Q, the dynamics of the 



III. THEORY 



Using the Cooper-pair resonant states representation 
one can recast the Hamiltonian JQl in the form 



H(t) = Ejoj, sin Qt + ea\ cos(ujt + ip) , 



(2) 



where e = 2eSQo/C and Ui are the Pauli matrices. As 
already anticipated, we will assume that the following 
relations arc fulfilled: 



e, Ml < huj < Ej < A , 

w/fi = N +p/q,p < q,N ^> 1. 



(3) 



Here N, p, and q are integer numbers. The final result, 
other than details of calculations, does not depend in any 
essential way on the arithmetic properties of N,p, q, and 
for definiteness we take q to be an even number (p is 
odd). Under these conditions the Hamiltonian J2J) is a 
periodic function of time with period T = qTo. 
The dc Josephson current is given by expression 



1 f T 2 f T/2 

I=Tj dtTrp(t)3(t) = j? J o dtTi 



mm , (4) 



where 



2e dH{t) _ eEj 



(73 COS VLt 



is the current operator, p(t) is the density matrix of the 
two-level system, while its time evolution is governed by 
the Liouville-von Neumann equation 



dp 



[H, p] - lhv{p - peg)- 



(5) 



Here p eq (t) = (1 — er 3 tanh(/3 Ej sin fit) )/2 is the qua- 
sistatic density matrix for the unperturbed Hamiltonian 
(e = 0), (3 is the inverse temperature, v is the relaxation 
rate. The formal solution of Eq. (JSJ can be expressed as 



p{t) = v I dt' e^'-^uit, t')p eq (t')uHt, t') 
where 



(6) 



Hh,h) = Texp 



dtH{t) 
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is an evolution operator possessing the symmetry prop- 
erties 

u{h + 9To/2, h + qTo/2) = a 3 u{t 2 ,h)a 3 , (7) 

which derives from the symmetry of the Hamiltonian J5J , 

H(t + T/2) = a 3 H{t)a 3 . 

It follows that the density matrix p(t) satisfies the rela- 
tion 



pit + T/2) = a 3P (t)a 3 



(8) 



which we used in Eq. 10} . In the limit vT <C 1 this 
relation together with Eq. © immediately results in the 
equation 



(9) 



for the density matrix p = p(t = 0). Here v = vT/2 <C 1 
and 



U = a 3 u(qT /2,0) 



(10) 



F=% I ' ' dt a 3 u{T/2,t)p eqv {t)u\T/2,t)a 3 .(ll) 
1 Jo 

Equation © and the spectral properties of the unitary 
operator U will play a central role as we proceed. 

In the limit v = vT/2 <C 1 the solution of Eq.@ has 
the form 

p = I/2 + B(|Ai)<Ai| - |A 2 )(A 2 |) +O(0), 

where |Afe=i, 2 ) are the eigenvectors of the operator U with 
eigenvalues exp(iAfc) and B = ((Ai |F|Ai) - (A 2 |F|A 2 ))/2 . 
From the symmetry of the Hamiltonian it follows that 
a 2 H(t)a 2 = -H(t), |A 2 ) = ct 2 |Ai) and Ai + A 2 = n. The 
expression HXk/T gives the quasienergies of the periodic 
Hamiltonian, H(t); we will refer to A& as a quasienergy 
factor. 

The average Josephson current, Eq.Q, can now be 
cast in the form 



i = Bj2(-V k+1 (*k\io\\k), 



(12) 



2eE r T / 2 

I = — -f / dtcosQtuUt,0)a 3 u(t,0) . (13) 
KT Jo 

In order to calculate the population coefhcient B and the 
average current in the quasienergy states, (A&| 2o|Afc), it is 
convenient to introduce an effective Hamiltonian, H (t; £), 
of the form 



H(t; £) = H(t) + £cr 3 tanh((3Ej sin Sit) . 



(14) 



One can show that Eq. l|12fl for the Josephson current 
can be expressed as 



2ehuj d\i d\i 
I = hm — — 

irqT i^o dip <9£ 



(15) 



where Ai (<£>;£) is the appropriate quasienergy factor of 
the effective Hamiltonian (|14l) . 

Only resonant inter-level transitions are of importance 
if the microwave interaction energy, e, is much smaller 
than the characteristic inter-level distance 2Ej ~ %uj. 
Ignoring the nonresonant part of the Hamiltonian H(t; £) 
(i.e. using the resonance approximation) one finds that 
the evolution operator u((n+l)To/2, nTg/2) has the form 

u((n+ l)T /2,nT /2;£) = a ^ e - i ^ n+1 ^ T ^ 2+ ^' 7 ^ 2 

X U res (T /2, 0; e i{™T /2+ v )a 3 /2 a n ^ 

Here w res (To/2, 0; £) is the evolution operator generated 
by the Hamiltonian iJ res (t;£), where 

H res (t; £) = (Ej sin Of - hu/2 + 

£tanh([3Ejsm.Slt))(T 3 + so\ . 

From the symmetry of the Hamiltonian -ff res (f;£) it fol- 
lows that the operator « res (To/2, 0; £) has the form 

u res (T /2,0;O = 
e ie<T3 (^l - r 2 + iT(j 1 )e ie ' 73 = e l9,T3 cr 1 Le w < J3 . 

The parameter t is the resulting probability amplitude 
for an inter-level transition when the superconducting 
phase difference changes from — 7r to 0. 

Equation 116f) permits us to represent the desired evo- 
lution operator U of Eq. I|10[l as 



±ie' 



9-1 

£r 3 IJ«((n + l)To/2,nT /2;0 

n=0 



(17) 



-ia 3 (e-ip/2) 



ia 3 (9- V /2) 



where a = p/q and the sign (±) depend on the arithmetic 
properties of the natural numbers iV,p,and q. 

Now one can see that the eigenvalues of the unitary op- 
erator U(£) can be found by solving the difference equa- 
tion 



|V(n + l)) = Le l( ™ n+ ^\^{n)) 



(18) 



with the boundary condition \ip(q)) = e lXk \ip(0)}, k — 
1, 2. A " quasiperiodic" equation of this type may be an- 
alyzed by method developed in Ref. The result for 
the quasienergies A& from such an analysis is 

Afc = {-l) k a,rccos(T q cos qtp) + c k {p,q,N) , (19) 

where the constants Ck are irrelevant. 

Taking into account Eqs. (|15fl and l|19|) the average 
current can now be rewritten as 



I =- 



zfiuj T 2 ( q ^ sin 2qip ( dr 



2-kTq 1 - r 2 ? cos 2 qip \d£ J e=0 



■ (20) 
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In this work we restrict ourselves to the low tempera- 
ture limit, where 0Ej 1. In this limit H res (t;^) — 
H r es{t;£, — hui/2) and, as a consequence, <9r 2 /<9£ = 
—2h~ 1 dT 2 /duj. Therefore, in order to calculate the cur- 
rent at low temperatures one can use the Hamiltonian 

As already mentioned, only resonant inter- level tran- 
sitions are important here. In the time interval (0,Xo/2) 
the resonant condition 2Ej sin flU — Huj is fulfilled only 
at the two times tq = il^ 1 &rcsinhuj/2Ej and t\ — 
To/ 2 — to. As these "points" in time are reached, inter- 
level Landau-Zener transitions^ (see Fig. © occur with 
probability amplitude u; 1 / 2 e J,s ', where 



w = 1 



exp 



ire 



( hui 



\2Ej 



Dividing time interval (0, T /2) into periods of free evo- 
lution (where interactions with the microwave field can 
be neglected) and describing the passage through the res- 
onant points by a Landau-Zener scattering matrix, one 
finds the probability for an inter-level transition to be 





! t 


♦ 




! * U ! 





10.0 10-5 11.0 11.5 12.0 12.5 13.0 13.5 

co/n 



FIG. 4: Microwave-induced current I in units of Jo = eio/ir 2 , 
plotted as a function the microwave frequency uo normal- 
ized to Q = ujj/2. The result was obtained with ip = 
7r/ll, hu/2Ej — 0.001, w = 0, 5 and q ma x = 7 (see text). 



IV. CONCLUSION 



r = 2w\/ 1 — w 2 cos(6 + d) 

pTo/2 — to 

Q = h~ 1 dtEjsmflt- uj(T /i-t ) 

J t 

and the final expression for the microwave- induced dc 
Josephson current becomes 

T eto ( hu\ t 2 i tan(6 + tf) . „ 

1 = arccos ^ ^ sm2qw . (21 

7T 2 V 2Ej J 1 - t 2 ? cos 2 qip r v ; 

Equation (|21|l determines the current-voltage character- 
istics of the SCPT studied, which depend crucially on 
the frequency of the microwave irradiation. Its validity 
is restricted to the special values of lo and V determined 
by Eq. © 

Averaging the current over a long time T (correspond- 
ing to a large value of g, since T = qTo), and plotting 
it as a function of the parameter w/VL, results in many 
sharp features. (Notice that two values of uj/U, which are 
arbitrarily close to each other might correspond to widely 
different values of q in Eq.©). It is clear, however, that 
any upper limit T max for the allowed averaging time will 
tend to smooth out the sharp features. The value of 
Tmax might be determined by either a relaxation time or 
a finite observation time. In either case currents calcu- 
lated for two values of (w /fi)i,2 will be indistinguishable 
if \q± — <?2 1 > QTmax- This argument provides a justifica- 
tion for smoothening the current as given by Eq. (|21|) in 
order to get a regular dependence / = I(uj/Q). An ap- 
propriate procedure is to evaluate Eq. lj2"T)l at all rational 
points for which q < q ma x- Then a smooth dependence 
can be obtained by making use of an interpo- 
lation procedure. The result of such an interpolation 
procedure is presented in Fig. 0] 



We have shown that the microwave dynamics of a 
single-Cooper-pair transistor is strongly influenced by in- 
terference effects caused by the quantum dynamics of 
the superconducting dot, which — as in Fig. [21 — forms 
the central island in the transistor structure. The sharp 
peaks that were found in the current-voltage (I — V) 
characteristics for "fractional" values of the voltage, as 
shown in Fig. 0] are a signature of a resonant interac- 
tion with the microwave field. The peaks occur when 
eV = huj/(N + p/q), where A, p, and q are integers; 
they are finite if N — > oo and do not vanish in the 
limit of weak microwave irradiation, The interaction with 
the microwave field results in multiple, coherent transi- 
tions between resonant Cooper-pair levels. This picture, 
which differs qualitatively from the Shapiro effecti, is a 
direct manifestation of the role of strong Coulomb cor- 
relations in the nonequilibrium superconducting dynam- 
ics of mesoscopic weak links. The Coulomb blockade of 
Cooper pair tunnelling reduces the number of relevant 
charge states to two, and is responsible for the ensuing 
quantum coherent two-level dynamics of Cooper pair res- 
onant states. 

A single-Cooper-pair transistor requires the Coulomb 
blockade energy E q to be much larger than the Josephson 
coupling energy Ej. If instead E j 3> E q , a large num- 
ber of dot charge states with different number of Cooper 
pairs are degenerate. As a result tunnel coupling to the 
leads form dot states with a well defined superconducting 
phase. Microwave properties in this limit are governed 
by the classical dynamics of the superconductiong phases, 
caused by the driving voltage and the microvawe field. If 
we were to imagine that E q could be diminished from a 
value much larger than Ej to a value much smaller than 
Ej, we would witness a transition from a situation where 
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the quantum dynamics of the superconducting phase was 
relevant to one where the classical dynamics would dom- 
inate. During such a transition the resonance microwave 
absorption would disappear and be replaced by standard 
Shapiro signatures, known from macroscopic weak links. 
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